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Deconfining transition in Full QCD 
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We present evidence that in full QCD with two dynamical quarks confinement is produced by dual supercon- 
ductivity of the vacuum as in the quenched theory. Preliminary information is obtained on the nature of the 
deconfining transition. 
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1. Introduction 

A schematic phase diagram of full QCD with 
2 dynamical flavours (m„ = = Wg) is shown 
in Fig. 1. The upper part of the diagram is the 
deconfined phase, the lower part is confined. The 
line is determined by the maxima of a number of 
susceptibilities ||l|,||, including the susceptibility 
XL of the Polyakov line 



XL 



d^a;(L(f,0)Lt(0,0)) 



(1) 



and the susceptibility x^h of the chiral order pa- 
rameter 



xch = / rf'x(v;^(f , o)v;V'(o, 0)) 



(2) 



All of them have a maximum at the same value 
of r, for a given m^, which defines the line in 
Fig. 1. For niq > 3 GeV, the maxima of xl di- 
verge proportionally to the volume V , indicating 
a first order transition. At rriq = there are theo- 
retical reasons and numerical indications that the 
transition is second order [pP,01 . At intermediate 



values of niq (tiny part of the line in Fig. 1) the 
susceptibilities do not diverge with V , and this is 
interpreted as absence of a phase transition: the 
line would correspond to a crossover. 
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Across the transition the density of the free en- 
ergy is a function of the order parameters. The 
singularities of derivatives of F are related to sus- 
ceptibilities of the order parameters. In our case 
L is the order parameter only at — oo, tpTp 
only at ruq — 0. A good order parameter in the 
whole range of to^'s could lead to a different as- 
signment for the order of the transition. 

Such a parameter could be the disorder pa- 
rameter (/i) which describes condensation of mag- 
netic charges [Q. That parameter has been con- 
structed and tested in quenched theory and 
its definition can be extended to full QCD 
In the spirit of Nc oo arguments one expects 
that the mechanism of confinement is the same in 
quenched and full QCD, quark loops being non 
leading in 1/A^c expansion. The operator fj, cre- 
ates a magnetic charge, as defined by some given 
abelian projection. (^) ^ signals dual super- 
conductivity. In the quenched theory the specific 
choice of the abelian projection proves to be im- 
material |0,D . For the details about the definition 
of fi we refer to |,|,|,|. For the quenched the- 
ory it is found that (/i) ^ for T < Tc and that 
(/i) = for T > Tc, ~ r-^ as T ^ T^, with 
r = 1 — T/Tc- A finite size scaling analysis of 
the infinite volume limit yields 5 = 0.20(3) for 
SU{2); d = 0.50(2) for SU{3). The critical in- 
dex 1/ has the values v — 0.62(1) for SU{2), in 
agreement with Ref. 0, ^ = 0.33(1) for SU{3), 
corresponding to a first order transition pd| . We 
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Figure 2. p as a function of Ng for T < Tc- 



repeat the same analysis for 2 staggered flavours, 
using the Wilson action for the pure gauge sector, 
and 12^ x 4, 16^ x 4,32^ x 4 lattices. The machines 
used are APEmillc crates. Part of the results are 
already published in Ref . . 

2. Numerical Results 

As usual instead of (/i) itself we determine the 
quantity P = log(/i) in terms of which {fi} = 

exp ( p{l3')dl3'^ . For T < Tc wc find that p is 
practically size independent, i.e. that (/i) ^ as 



Figure 3. p as a function of Ns for T > Tc- As 
A'^s — > oo, p — + — oo if the magnetic charge is non- 
zero and stays constant otherwise. 



■ oo (see Fig. 2). For T > we find (Fig. 3) 



p = —kNg + const. 



(3) 



i.e. that (/i) is strictly zero as A^s — > oo. 
More extended checks of superselection of mag- 
netic charge in the deconfined phase have been 
done, showing that for different magnetic charges 
created by p, p ^ — oo if the net magnetic charge 
is non-zero, p — > const, if it is zero. Some exam- 
ples are given in Fig. 3. 

Around Tc a finite size scaling analysis goes 
as follows. By dimensional arguments one can 
parametrize (p) as 



(4) 



where a is the lattice spacing and t = I — T/Tc- 
If the correlation length goes large, <^ ~ t^"^, 
a/^ 1 and the dependence on a can be ne- 
glected. The variable A^s/C can be traded for 
tNs^'^ ■ In the quenched theory niq is absent and 
(/i) = t^4>{0,tNI^''), which gives for p the scal- 
ing behaviour p/Ns^" = /{tN}^"), that can be 
tested and allows the determination of T^, v and 
6. In particular it implies that the height of the 



peak scales as A, 
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Figure 4. p peak at different values of and 
TOTr/mp ~ 0.505 

In presence of dynamical quarks, vriq sets 
another scale, and we face a two scale prob- 
lem. Simulations done at fixed {ttlt^ /nip) |^ show 
that the height of the peak p roughly scales as 
= V (see Fig. 4). A more refined analysis 
can be done by choosing ruq and Ng such that 
niqN] in Eq. (^) is kept constant, and the prob- 
lem is again reduced to a one scale problem, al- 
lowing to determine v, which gives information 
on the order of the transition. The index 7 is 
known to be 7 ~ 2.49 so that the three sets 



{mq = 0.075, iV^ = 16), (m. 



0.043, Ns 



20) 



and {niq = 0.01335, Ng = 32) keep mqN] con- 
stant. Preliminary results indicate that v = 1/3, 
or that the transition is compatible with first or- 
der. This is clearly visible from Fig. 5, where the 
height of the peak is plotted versus N^. More de- 
tailed numerical analyses involving more values 



oiNg 



and the use of an improved action, to 



improve the quality of scaling, are needed to draw 
a definite conclusion. 

In any case it is demonstrated that also in full 
QCD confinement is produced by monopole con- 
densation. 
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Figure 5. Values of p at the peak as a function 
of nI'"", with V = l/'i. The height of the peak 
clearly scales as A''^''. 



